The solar oxygen crisis: Probably not the last word 



O 

o 

(N 



H. Socas-Navarro 

High Altitude Observatory, NCAR] 3080 Center Green Dr, Boulder, CO 80301, USA 

navarroOucar . edu 

A. A. Norton 

National Solar Observatory, 950 N. Cherry Ave., Tucson, AZ 85719, USA 

nortonOnso . edu 



(N 

> 
(N 

(N 
O 
l> 

o 
o 



ABSTRACT 

In this work we present support for recent claims that advocate a downward revision of the 
solar oxygen abundance. Our analysis employs spatially-resolved spectro-polarimetric observa- 
tions including the Fe I lines at 6302 A and the O I infrared triplet around 7774 A in the quiet 
Sun. We used the Fe I lines to reconstruct the three-dimensional thermal and magnetic structure 
of the atmosphere. The simultaneous O I observations were then employed to determine the 
abundance of oxygen at each pixel, using both LTE and non-LTE (NLTE) approaches to the 
radiative transfer. In this manner, we obtain values of logeo=8.63 (NLTE) and 8.93 (LTE) dex. 
We find an unsettling fluctuation of the oxygen abundance over the field of view. This is likely 
an artifact indicating that, even with this relatively refined strategy, important physical ingredi- 
ents are still missing in the picture. By examining the spatial distribution of the abundance, we 
estimate realistic confidence limits of approximately 0.1 dex. 

Subject headings: Sun: abundances - Sun: granulation - stars: abundances - stars: atmospheres 



Oxygen is the third most abundant chemical el- 
ement in the universe. Its abundance in the solar 
atmosphere has been traditio nally accepted to be 
appro ximately logeo=8.93 (jAnders fc Grevesse 
with only minor co rrections until re- 
centl y (e.g., logeo=8.83 by 'Grevess e fc Sauval 
19981 ). These determinations are based on the 
use of one-dimensiona l (ID) semiempirical mod- 
els (such as those o f 'Holweger & Mueller' I l974 ; 



20061: IScott et al.ll2006l ). 



2004 : iMelended 120041) and CO (e.g.. lAvres et al 



Co ntroversy was sparked when lAsplund et al 



Gingerich et al.lll971t tVernazza et al. 1981) to fit 



spectral features produced by oxygen atoms and 
molecules. The most common abundance indi- 
cators are: the forb i dden [O I] line at 6300 A 
(e.g., iLambertl Il978t lAllende Prieto et al.ll200lh . 
the infrared ( IR) t ri plet around 7774 A (e.g., 



lAsplund et all 120041: Ishchukina et all l2005l) 



well as molecular bands of OH (e.g. . lAsplund et al 



as 



(j2004[ ) applied a new hydrodynamic simulation 
of the solar granulation to the determination of 
the oxygen abundance, resulting in a consider- 
ably lower valued of logeo=8.66. The new value 
presents some advantages, such as better agree- 
ment between various indicators (IR triplet, [O I] 
and OH vibrational and roto- vibrational bands). 
It also leads to a better fit of the solar chemical 
com position within its ga lactic environment (see, 
e.g., ISofia fc Meveiil200ll and references therein). 
On the down side, the traditional composition 
had led to a truly remarkable agreement between 
the sound speed predicted by solar interior mod- 



^The National Center for Atmospheric Research 
(NCAR) is sponsored by the National Science Foundation. 



^ A lthough a low abund ance had been previously derived 
on lAUende Prieto et al.| 1200 li ) based on the analysis of the 
forbidden line at 6300 A. 
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els and the helioseismical inversions. A revision 
of the oxygen abundance, along with the cascad- 
ing effects on nitrogen, carbon and neon, would 
require a thorough rew orking of that comparison 
and some recent work ( Delahave &: Pinsonneaul'3 
l2006t iBasu et aL 2007 ) already argues for its in- 
compatibility. The far-reaching implications of 
the proposed change (notice that it represents a 
factor of two in the actual n umber densities) has 
prompted Avres et al.l ( 20061 ) to label this problem 
as the solar oxygen crisis. 

(|2006h derive a 



sity matches that of the Harvard S mithsonian Ref- 



erenc e Atmosphere (HSRA, see iGingerich et al 



19711 ). Although microturbulence is permitted as 



In their work, lAvres et al 



semiempirical ID model of the average thermal 
stratification and use it to obtain the oxygen abun- 
dance from observed CO bands. This approach 
led them to recommend a rather high value of 
logeo=8.85, consistent with the traditional de- 
terminati ons. While di s cussin g the discrepancy 
with the lAsplund et al.l (|2004l ) results, they cor- 
rectly point out that a 3D theoretical simulation 
is not necessarily superior to a ID semiempirical 
model when used for the diagnostics of observa- 
tions. They both have advantages and drawbacks. 

The approach taken here tries to capture the 
best elements of both worlds, by building a 
semiempirical 3D model of the solar photosphere. 
The data set was acquired on 2006 Oct 27 with 
the Spectro-Polarime ter for Infrared and Opt i- 
cal Regions (SPINOR. ISocas-Navarro et al.l l2006V 
It consists of a high-resolution (0.7") spectro- 
polarimetric scan of a quiet Sun region including 
two Fe I lines at 6302 A and the oxygen infrared 
triplet at 7774 A. The Fe I lines have been used to 
build the semiempirical 3D model, while the O I 
triplet is employed for the abundance determina- 
tion. Since both sets of lines have been observed 
simultaneously, there is full consistency between 
the data employed to construct the model and to 
derive the abundances. 

Our 3D model has been obtained from the in- 
version of the well-known Fe I lines at 6301.5 (for 
which the oscillat or strength has b een measured 
in the laboratorv. iBard et al. 1991 ) and 6302.5 A 
together with the continuum at 7774 A. For each 
pixel in the field of view, we derived the vertical 
stratification of temperature, density, line-of-sight 
velocity and longitudinal magnetic field (details on 
the Fe I observations and the inversion will be pro- 
vided elsewhere). We impose the condition that 
the average observed quiet Sun continuum inten- 



a free parameter, the inversion code always re- 
trieved very small values (~100 m s~^). The 
model is publicly available upon request from the 
authors. The spatial resolution is obviously lim- 
ited when compared to that of the Asplund et al. 
simulations, but it is sufficient to resolve the gran- 
ulation. Moreover, those pixels exhibiting sigin- 
ficant polarization signal have been treated with 
a two-component scenario (one for the magnetic 
concentration and another for the non-magnetic 
surroundings) . Fig [T] shows a composite image of 
the temperature and the magnetic flux density at 
the base of the photosphere. 

The oxygen abundance was determined inde- 
pendently at each pixel of the map, using the 
model atmosphere from the Fe I inversion, and ac- 
counting for magnetic fields (i.e., considering the 
Zeeman splitting and solving the Stokes vector 
transfer equation) where appropriate. Synthetic 
profiles were computed at intervals of 0.1 dex, 
and the x^"VS-logeo curve was interpolated to 
find the minimum with '^O.Ol dex accuracy (the 
curves obtained are very smooth). For the 
LTE calcul ations we used the Stokes synthesis 
code LILIA (|Socas-Navarroll200l[ ). whereas for the 
non-LTE case we employed the code developed by 
Socas-Navarro et al. (2000). To alleviate the com- 
putational burden we used a 1.5D approach (i.e., 
neglecting lateral radiative tra nsfer). This approx- 
imation is justified by Asplun d et al.l ((2004), who 
obtained very similar results when comparing 3D 
and 1.5D syntheses in their simulation. The model 
atom considered has 13 O I levels plus the con- 
tinuum, with 18 bound-bound transitions and 13 
bound-free (6 of them treated with radiation tem- 
peratures). The atomic data are the same as in 
Carlsson fc Judge! (|l993h . 



Fo llowing the recommendation of lAsplund et al. 
(|2004) , we neglect excitation and ionization due to 
inelastic collisions with neutral hydrogen. For the 
collisional broadenin g of t he lin es we employed 
the formulation of Gravl ( 1976f ) and tested the 
influence of the Van der Waals enhancement pa- 
rameter on the oxygen abundance. For this test 
we took the spatially-un resolved disk center atlas 
of iNeckel fc Labs! (| 19841 ) . inverted the Fe I pair 
of lines around 6302 A and used the resulting ID 
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model to synthesize the O I IR triplet. Vary- 
ing the damping enhancement parameter between 
and 10 had very little influence on the oxygen 
abundance (~0.01 dex). For larger values, the syn- 
thetic profiles are too broad to fit the observations. 
While the enhancement parameter strongly influ- 
ences the equivalent widths of the lines (because it 
broadens the wings leaving the cores more or less 
intact), its effect is less pronounced in our analysis 
that considers the least-squares fit to the overall 
profile. In this test with unresolved observations 
and a ID model we obtained an oxygen abundance 
of log eo=8.66, which is consistent with the pre- 



vious ID NLTE determinat i ons of lAsplund et al 
(|2004l ) and lShchukina et al.l (|2005l ). 



Since the O I lines are shallower than the Fe I 
lines, one would intuitively expect their formation 
height range to be narrower. In other words, the 
height range that we need for the synthesis of the 
O I triplet is contained entirely in the range that 
we can reconstruct from the Fe I lines. A formal 
verification of this is presented in Fig [2l which 
sho ws the tem peratu re response functions (RFs, 



see 



iandi Deg PInno centi fc Landi Degl'Innocenti 



19771 ) of the lines at several wavelengths, from 
the continuum to the line center. RFs measure 
how the line profile reacts to atmospheric pertur- 
bations at various heights, and therefore provide 
very valuable information on the sensitivity range 
of spectral features. 

The spatial distribution of log eo obtained with 
our 3D semiempirical model is shown in Fig[3l In 
agreement with previous works, we find that the 
LTE approximation largely overestimates logeo 
by approximately 0.3 dex (this value is slightly de- 
pendent on spatial location). The variation of the 
inferred log eo over the field of view is due to fluc- 
tuations in the solar atmospheric conditions. We 
can see that the abundance in the pore is consid- 
erably higher, reaching values up to 0.3 dex larger 
than the average, and decreasing gradually as we 
move away from it. In general, magnetic concen- 
trations exhibit higher abundances. This suggests 
that previous abundance determinations might be 
slightly biased by solar (or, more generally, stellar) 
activity, depending on the phase of the cycle in 
which the observations were taken. The mean val- 
ues that we obtained over the entire fleld of view 
are logeo=8.94 (LTE) and logeo=8.64 (NLTE) 
dex, with standard deviations of 0.08 dex. In the 



NLTE case, this is equivalent to a concentration 
of eo/e// of 442 (±80) xlQ-^. If we restrict the 
analysis to pixels with less than 100 G of flux den- 
sity, we obtain logeo=8.93 (LTE) and logeo=8.63 
(NLTE). 

Obviously, there is no physical reason to ex- 
pect the actual abundance to exhibit spatial vari- 
ations in the solar photosphere. We must then 
conclude that this is an artifact of the analysis, 
probably due to imperfect modeling especially in 
the presence of magnetic fields (notice that the 
granulation pattern is not visible in the abundance 
images). We stress that the modeling employed 
here should be more reliable than that in previ- 
ous works (since it considers 3D geometry, mag- 
netic fields, NLTE and is based on actual obser- 
vations). Yet, it still falls short of providing a 
solid, accurate, abundance determination. The 
uncertai nties (^0.1 dex) are larger than pr eviously 
thought (*Pinson neault &: Delahav^ 20061 have al- 
ready pointed this out and their error budget is 
consistent with our figure), and this is also an im- 
portant point to bear in mind. 
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Fig. 1. — Maps of the temperature and the mag- 
netic flux density at the base of the photosphere in 
the 3D model, as derived from the SPINOR data. 
Left panel: Temperature. In those pixels that 
exhibit polarization signal, the average tempera- 
ture inside and outside the magnetic concentration 
(weighted according to occupation fraction) is de- 
picted. The cool feature on the lower left corner 
is a pore that was used to stabilize the adaptive 
optics, allowing us to achieve the spatial resolu- 
tion shown in the image. Right panel: Decimal 
logarithm of the magnetic flux density in G. 





Fig. 2. — Normalized response functions to tem- 
perature perturbations, computed in the HSRA 
quiet Sun model. Left panels: Line profiles of the 
Fe I 6301 A (top) and O I 7772 A (bottom) lines. 
Right panels: Response of the intensity profile at 
the wavelengths marked with vertical dotted lines 
in the left panels, as a function of height where 
the perturbation occurs. The 7772 A line is the 
deepest (darkest core) of the multiplet and there- 
fore spans the broadest height range of the O I 
triplet. Nevertheless, its formation range is still 
well within that spanned by the Fe I lines. 
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Fig. 3. — Logarithmic oxygen abundance logeo 
(with the usual convention in Astrophysics that 
logeH=12). Left panel: LTE calculation. Right 
panel: NLTE calculation. Note that the color 
scales are different, since the LTE abundances are 
systematically higher. The vertical lines are miss- 
ing data due to errors in the detector system that 
was recording the oxygen data. 



